SplicerAV: a tool for mining microarray expression data for changes in RNA processing.

BACKGROUND: Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets. RESULTS: Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets. CONCLUSIONS: Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

Traumatic brain injury exacerbates neurodegenerative pathology: improvement with an apolipoprotein E-based therapeutic.

Cognitive impairment is common following traumatic brain injury (TBI), and neuroinflammatory mechanisms may predispose to the development of neurodegenerative disease. Apolipoprotein E (apoE) polymorphisms modify neuroinflammatory responses, and influence both outcome from acute brain injury and the risk of developing neurodegenerative disease. We demonstrate that TBI accelerates neurodegenerative pathology in double-transgenic animals expressing the common human apoE alleles and mutated amyloid precursor protein, and that pathology is exacerbated in the presence of the apoE4 allele. The administration of an apoE-mimetic peptide markedly reduced the development of neurodegenerative pathology in mice homozygous for apoE3 as well as apoE3/E4 heterozygotes. These results demonstrate that TBI accelerates the cardinal neuropathological features of neurodegenerative disease, and establishes the potential for apoE mimetic therapies in reducing pathology associated with neurodegeneration.

Microbial inactivation of Pseudomonas putida and Pichia pastoris using gene silencing.

Antisense deoxyoligonucleotide (ASO) gene silencing was investigated as a potential disinfection tool for industrial and drinking water treatment application. ASOs bind with their reverse complementary mRNA transcripts thereby blocking protein translation. While ASO silencing has mainly been studied in medicine, it may be useful for modulating gene expression and inactivating microorganisms in environmental applications. In this proof of concept work, gene targets were sh ble (zeocin resistance) and todE (catechol-2,3-dioxygenase) in Pichia pastoris and npt (kanamycin resistance) in Pseudomonas putida. A maximum 0.5-fold decrease in P. pastoris cell numbers was obtained following a 120 min incubation with single-stranded DNA (ssDNA) concentrations ranging from 0.2 to 200 nM as compared to the no ssDNA control. In P. putida, a maximum 5.2-fold decrease was obtained after 90 min with 400 nM ssDNA. While the silencing efficiencies varied for the 25 targets tested, these results suggest that protein activity as well as microbial growth can be altered using ASO gene silencing-based tools. If successful, this technology has the potential to eliminate some of the environmental and health issues associated with the use of strong chemical biocides. However, prior to its dissemination, more research is needed to increase silencing efficiency and develop effective delivery methods.

MRP3: a molecular target for human glioblastoma multiforme immunotherapy.

BACKGROUND: Glioblastoma multiforme (GBM) is refractory to conventional therapies. To overcome the problem of heterogeneity, more brain tumor markers are required for prognosis and targeted therapy. We have identified and validated a promising molecular therapeutic target that is expressed by GBM: human multidrug-resistance protein 3 (MRP3). METHODS: We investigated MRP3 by genetic and immunohistochemical (IHC) analysis of human gliomas to determine the incidence, distribution, and localization of MRP3 antigens in GBM and their potential correlation with survival. To determine MRP3 mRNA transcript and protein expression levels, we performed quantitative RT-PCR, raising MRP3-specific antibodies, and IHC analysis with biopsies of newly diagnosed GBM patients. We used univariate and multivariate analyses to assess the correlation of RNA expression and IHC of MRP3 with patient survival, with and without adjustment for age, extent of resection, and KPS. RESULTS: Real-time PCR results from 67 GBM biopsies indicated that 59/67 (88%) samples highly expressed MRP3 mRNA transcripts, in contrast with minimal expression in normal brain samples. Rabbit polyvalent and murine monoclonal antibodies generated against an extracellular span of MRP3 protein demonstrated reactivity with defined MRP3-expressing cell lines and GBM patient biopsies by Western blotting and FACS analyses, the latter establishing cell surface MRP3 protein expression. IHC evaluation of 46 GBM biopsy samples with anti-MRP3 IgG revealed MRP3 in a primarily membranous and cytoplasmic pattern in 42 (91%) of the 46 samples. Relative RNA expression was a strong predictor of survival for newly diagnosed GBM patients. Hazard of death for GBM patients with high levels of MRP3 RNA expression was 2.71 (95% CI: 1.54-4.80) times that of patients with low/moderate levels (p = 0.002). CONCLUSIONS: Human GBMs overexpress MRP3 at both mRNA and protein levels, and elevated MRP3 mRNA levels in GBM biopsy samples correlated with a higher risk of death. These data suggest that the tumor-associated antigen MRP3 has potential use for prognosis and as a target for malignant glioma immunotherapy.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

An integrated approach to the prediction of chemotherapeutic response in patients with breast cancer.

BACKGROUND: A major challenge in oncology is the selection of the most effective chemotherapeutic agents for individual patients, while the administration of ineffective chemotherapy increases mortality and decreases quality of life in cancer patients. This emphasizes the need to evaluate every patient's probability of responding to each chemotherapeutic agent and limiting the agents used to those most likely to be effective. METHODS AND RESULTS: Using gene expression data on the NCI-60 and corresponding drug sensitivity, mRNA and microRNA profiles were developed representing sensitivity to individual chemotherapeutic agents. The mRNA signatures were tested in an independent cohort of 133 breast cancer patients treated with the TFAC (paclitaxel, 5-fluorouracil, adriamycin, and cyclophosphamide) chemotherapy regimen. To further dissect the biology of resistance, we applied signatures of oncogenic pathway activation and performed hierarchical clustering. We then used mRNA signatures of chemotherapy sensitivity to identify alternative therapeutics for patients resistant to TFAC. Profiles from mRNA and microRNA expression data represent distinct biologic mechanisms of resistance to common cytotoxic agents. The individual mRNA signatures were validated in an independent dataset of breast tumors (P = 0.002, NPV = 82%). When the accuracy of the signatures was analyzed based on molecular variables, the predictive ability was found to be greater in basal-like than non basal-like patients (P = 0.03 and P = 0.06). Samples from patients with co-activated Myc and E2F represented the cohort with the lowest percentage (8%) of responders. Using mRNA signatures of sensitivity to other cytotoxic agents, we predict that TFAC non-responders are more likely to be sensitive to docetaxel (P = 0.04), representing a viable alternative therapy. CONCLUSIONS: Our results suggest that the optimal strategy for chemotherapy sensitivity prediction integrates molecular variables such as ER and HER2 status with corresponding microRNA and mRNA expression profiles. Importantly, we also present evidence to support the concept that analysis of molecular variables can present a rational strategy to identifying alternative therapeutic opportunities.

Gene expression disruptions of organism versus organ in Drosophila species hybrids.

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

A switch in the control of growth of the wing imaginal disks of Manduca sexta.

BACKGROUND: Insulin and ecdysone are the key extrinsic regulators of growth for the wing imaginal disks of insects. In vitro tissue culture studies have shown that these two growth regulators act synergistically: either factor alone stimulates only limited growth, but together they stimulate disks to grow at a rate identical to that observed in situ. It is generally thought that insulin signaling links growth to nutrition, and that starvation stops growth because it inhibits insulin secretion. At the end of larval life feeding stops but the disks continue to grow, so at that time disk growth has become uncoupled from nutrition. We sought to determine at exactly what point in development this uncoupling occurs. METHODOLOGY: Growth and cell proliferation in the wing imaginal disks and hemolymph carbohydrate concentrations were measured at various stages in the last larval instar under experimental conditions of starvation, ligation, rescue, and hormone treatment. PRINCIPAL FINDINGS: Here we show that in the last larval instar of M. sexta, the uncoupling of nutrition and growth occurs as the larva passes the critical weight. Before this time, starvation causes a decline in hemolymph glucose and trehalose and a cessation of wing imaginal disks growth, which can be rescued by injections of trehalose. After the critical weight the trehalose response to starvation disappears, and the expression of insulin becomes decoupled from nutrition. After the critical weight the wing disks loose their sensitivity to repression by juvenile hormone, and factors from the abdomen, but not the brain, are required to drive continued growth. CONCLUSIONS: During the last larval instar imaginal disk growth becomes decoupled from somatic growth at the time that the endocrine events of metamorphosis are initiated. These regulatory changes ensure that disk growth continues uninterrupted when the nutritive and endocrine signals undergo the drastic changes associated with metamorphosis.

The effects of sleep hypoxia on coagulant factors and hepatic inflammation in emphysematous rats.

OBJECTIVES: To develop a sleep hypoxia (SH) in emphysema (SHE) rat model and to explore whether SHE results in more severe hepatic inflammation than emphysema alone and whether the inflammation changes levels of coagulant/anticoagulant factors synthesized in the liver. METHODS: Seventy-five rats were put into 5 groups: SH control (SHCtrl), treated with sham smoke exposure (16 weeks) and SH exposure (12.5% O(2), 3 h/d, latter 8 weeks); emphysema control (ECtrl), smoke exposure and sham SH exposure (21% O(2)); short SHE (SHEShort), smoke exposure and short SH exposure (1.5 h/d); mild SHE (SHEMild), smoke exposure and mild SH exposure (15% O(2)); standard SHE (SHEStand), smoke exposure and SH exposure. Therefore, ECtrl, SHEShort, SHEMild and SHEStand group were among emphysematous groups. Arterial blood gas (ABG) data was obtained during preliminary tests. After exposure, hepatic inflammation (interleukin -6 [IL-6] mRNA and protein, tumor necrosis factor α [TNFα] mRNA and protein) and liver coagulant/anticoagulant factors (antithrombin [AT], fibrinogen [FIB] and Factor VIII [F VIII]) were evaluated. SPSS 11.5 software was used for statistical analysis. RESULTS: Characteristics of emphysema were obvious in emphysematous groups and ABGs reached SH criteria on hypoxia exposure. Hepatic inflammation parameters and coagulant factors are the lowest in SHCtrl and the highest in SHEStand while AT is the highest in SHCtrl and the lowest in SHEStand. Inflammatory cytokines of liver correlate well with coagulant factors positively and with AT negatively. CONCLUSIONS: When SH is combined with emphysema, hepatic inflammation and coagulability enhance each other synergistically and produce a more significant liver-derivative inflammatory and prothrombotic status.

Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals.

There is great interindividual variability in HIV-1 viral setpoint after seroconversion, some of which is known to be due to genetic differences among infected individuals. Here, our focus is on determining, genome-wide, the contribution of variable gene expression to viral control, and to relate it to genomic DNA polymorphism. RNA was extracted from purified CD4+ T-cells from 137 HIV-1 seroconverters, 16 elite controllers, and 3 healthy blood donors. Expression levels of more than 48,000 mRNA transcripts were assessed by the Human-6 v3 Expression BeadChips (Illumina). Genome-wide SNP data was generated from genomic DNA using the HumanHap550 Genotyping BeadChip (Illumina). We observed two distinct profiles with 260 genes differentially expressed depending on HIV-1 viral load. There was significant upregulation of expression of interferon stimulated genes with increasing viral load, including genes of the intrinsic antiretroviral defense. Upon successful antiretroviral treatment, the transcriptome profile of previously viremic individuals reverted to a pattern comparable to that of elite controllers and of uninfected individuals. Genome-wide evaluation of cis-acting SNPs identified genetic variants modulating expression of 190 genes. Those were compared to the genes whose expression was found associated with viral load: expression of one interferon stimulated gene, OAS1, was found to be regulated by a SNP (rs3177979, p = 4.9E-12); however, we could not detect an independent association of the SNP with viral setpoint. Thus, this study represents an attempt to integrate genome-wide SNP signals with genome-wide expression profiles in the search for biological correlates of HIV-1 control. It underscores the paradox of the association between increasing levels of viral load and greater expression of antiviral defense pathways. It also shows that elite controllers do not have a fully distinctive mRNA expression pattern in CD4+ T cells. Overall, changes in global RNA expression reflect responses to viral replication rather than a mechanism that might explain viral control.

Cryptococcal cell morphology affects host cell interactions and pathogenicity.

Cryptococcus neoformans is a common life-threatening human fungal pathogen. The size of cryptococcal cells is typically 5 to 10 microm. Cell enlargement was observed in vivo, producing cells up to 100 microm. These morphological changes in cell size affected pathogenicity via reducing phagocytosis by host mononuclear cells, increasing resistance to oxidative and nitrosative stress, and correlated with reduced penetration of the central nervous system. Cell enlargement was stimulated by coinfection with strains of opposite mating type, and ste3aDelta pheromone receptor mutant strains had reduced cell enlargement. Finally, analysis of DNA content in this novel cell type revealed that these enlarged cells were polyploid, uninucleate, and produced daughter cells in vivo. These results describe a novel mechanism by which C. neoformans evades host phagocytosis to allow survival of a subset of the population at early stages of infection. Thus, morphological changes play unique and specialized roles during infection.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Impact of simian immunodeficiency virus infection on chimpanzee population dynamics.

Like human immunodeficiency virus type 1 (HIV-1), simian immunodeficiency virus of chimpanzees (SIVcpz) can cause CD4+ T cell loss and premature death. Here, we used molecular surveillance tools and mathematical modeling to estimate the impact of SIVcpz infection on chimpanzee population dynamics. Habituated (Mitumba and Kasekela) and non-habituated (Kalande) chimpanzees were studied in Gombe National Park, Tanzania. Ape population sizes were determined from demographic records (Mitumba and Kasekela) or individual sightings and genotyping (Kalande), while SIVcpz prevalence rates were monitored using non-invasive methods. Between 2002-2009, the Mitumba and Kasekela communities experienced mean annual growth rates of 1.9% and 2.4%, respectively, while Kalande chimpanzees suffered a significant decline, with a mean growth rate of -6.5% to -7.4%, depending on population estimates. A rapid decline in Kalande was first noted in the 1990s and originally attributed to poaching and reduced food sources. However, between 2002-2009, we found a mean SIVcpz prevalence in Kalande of 46.1%, which was almost four times higher than the prevalence in Mitumba (12.7%) and Kasekela (12.1%). To explore whether SIVcpz contributed to the Kalande decline, we used empirically determined SIVcpz transmission probabilities as well as chimpanzee mortality, mating and migration data to model the effect of viral pathogenicity on chimpanzee population growth. Deterministic calculations indicated that a prevalence of greater than 3.4% would result in negative growth and eventual population extinction, even using conservative mortality estimates. However, stochastic models revealed that in representative populations, SIVcpz, and not its host species, frequently went extinct. High SIVcpz transmission probability and excess mortality reduced population persistence, while intercommunity migration often rescued infected communities, even when immigrating females had a chance of being SIVcpz infected. Together, these results suggest that the decline of the Kalande community was caused, at least in part, by high levels of SIVcpz infection. However, population extinction is not an inevitable consequence of SIVcpz infection, but depends on additional variables, such as migration, that promote survival. These findings are consistent with the uneven distribution of SIVcpz throughout central Africa and explain how chimpanzees in Gombe and elsewhere can be at equipoise with this pathogen.

Hybrid transgenic mice reveal in vivo specificity of G protein-coupled receptor kinases in the heart.

G protein-coupled receptor kinases (GRKs) phosphorylate activated G protein-coupled receptors, including alpha(1B)-adrenergic receptors (ARs), resulting in desensitization. In vivo analysis of GRK substrate selectivity has been limited. Therefore, we generated hybrid transgenic mice with myocardium-targeted overexpression of 1 of 3 GRKs expressed in the heart (GRK2 [commonly known as the beta-AR kinase 1], GRK3, or GRK5) with concomitant cardiac expression of a constitutively activated mutant (CAM) or wild-type alpha(1B)AR. Transgenic mice with cardiac CAMalpha(1B)AR overexpression had enhanced myocardial alpha(1)AR signaling and elevated heart-to-body weight ratios with ventricular atrial natriuretic factor expression denoting myocardial hypertrophy. Transgenic mouse hearts overexpressing only GRK2, GRK3, or GRK5 had no hypertrophy. In hybrid transgenic mice, enhanced in vivo signaling through CAMalpha(1B)ARs, as measured by myocardial diacylglycerol content, was attenuated by concomitant overexpression of GRK3 but not GRK2 or GRK5. CAMalpha(1B)AR-induced hypertrophy and ventricular atrial natriuretic factor expression were significantly attenuated with either concurrent GRK3 or GRK5 overexpression. Similar GRK selectivity was seen in hybrid transgenic mice with wild-type alpha(1B)AR overexpression concurrently with a GRK. GRK2 overexpression was without effect on any in vivo CAM or wild-type alpha(1B)AR cardiac phenotype, which is in contrast to previously reported in vitro findings. Furthermore, endogenous myocardial alpha(1)AR mitogen-activated protein kinase signaling in single-GRK transgenic mice also exhibited selectivity, as GRK3 and GRK5 desensitized in vivo alpha(1)AR mitogen-activated protein kinase responses that were unaffected by GRK2 overexpression. Thus, these results demonstrate that GRKs differentially interact with alpha(1B)ARs in vivo such that GRK3 desensitizes all alpha(1B)AR signaling, whereas GRK5 has partial effects and, most interestingly, GRK2 has no effect on in vivo alpha(1B)AR signaling in the heart.